Measuring the Complete Transverse Spatial Mode Spectrum of a Wave Field 
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We put forward a method that allows the experimental determination of the entire spatial mode 
spectrum of any arbitrary monochromatic wave field in a plane normal to its propagation direction. 
For coherent optical fields, our spatial spectrum analyzer can be implemented with a small number 
of benchmark refractive elements embedded in a single Mach-Zehnder interferometer. We detail 
an efficient setup for measuring in the Hermite- Gaussian mode basis. Our scheme should also be 
feasible in the context of atom optics for analyzing the spatial profiles of macroscopic matter waves. 

PACS numbers: 42.25.-p, 42.30.-d, 42.15.Eq, 39.20.+q 



The temporal [1 , vectorial [2 and spatial [3 degrees 
of freedom of electromagnetic fields with increasing com- 
plexity are exploited in applications ranging from com- 
munications to medicine. The detailed knowledge of the 
spatial mode spectrum of these beams, generally emitted 
by laser devices, is fundamental in optimizing applica- 
tions as well as in research, both in classical and quan- 
tum optics. Many imaging and tailoring techniques, for 
instance, are based on the possibility to devise transfor- 
mations induced by optical elements by determining the 
effects on the spatial components of generic signals [4]. 
Recently, an intense research activity on multimode light 
has burgeoned in quantum information, communication 
and imaging [5 , with successful demonstrations span- 
ning from nanodisplacement measurements [6 , to par- 
allel information [7 , high-dimensional entanglement [8 , 
coherent transfer of suitable prepared superpositions of 
vortices onto Bose-Einstein condensates [9 , and hot Ru- 
bidium vapors [10 . A key ingredient in all these sce- 
narios is spatial encoding^ which conveys several inde- 
pendent channels of information by exploiting the access 
to a large number of optical transverse modes. The use 
of multimode beams thus demands for the development 
of techniques allowing to characterize their spatial spec- 
trum and to retrieve the information encoded in differ- 
ent components. Fourier modes are of immediate access 
through a basic lens setup [4 , while rotating elements 
lead to 'helical' spectra decomposition [11 . To the best 
of our knowledge, however, no direct method is known to 
measure the Hermite-Gaussian (HG) modes spectrum. 
Building on the symplectic formalism to describe first- 
order optical transformations, we present in this Letter 
a general strategy to find an arrangement of refractive 
elements that enables the quantitative measurement of 
the transverse spatial spectrum of light beams by using 
a single Mach-Zehnder interferometer. Due to the gener- 
ality of our approach, different experimental setups can 
be implemented to extract the full spatial spectrum of 



multimode transverse beams in many bases, including 
Laguerre- Gaussian (LG) modes. Here we focus on the 
spectrum analyzer for HG modes, as these are the most 
common in laser physics and appear naturally in devices 
where astigmatism, strain or slight misalignment drive 
the system toward rectangular symmetry [l2[. Further- 
more, the presented framework is rather suggestive of 
analog measurements for matter waves spatial spectra. 

Spatial modes are ubiquitous. Sets of modes u satis- 
fying the wave equation {idr^ + 9^ + = describe 
broad classes of physical systems. When the evolution 
variable is = 2/cz, the sets comprise the optical parax- 
ial modes [12 with wave number k propagating along 
the z direction. For r] = 2mt/h^ the sets model instead 
the quantum dynamics of free particles of mass m in the 
xy plane. Since there are two transverse spatial vari- 
ables, each entire set of mode solutions Um,n is labeled 
by two integers m, n. Relevant examples are the HG and 
LG mode bases. Any scalar field obeying the above 
wave equation can thus be decomposed in terms of these 
modes as = ^ Cm,nUm,n- A fundamental question 
then arises: How can one measure the mode spectrum 
(weights) Vm,n = |C'm,nP of any given scalar field? 

To answer the above question we apply a symplectic 
invariant approach [13l HH |T5] . Symplectic methods have 
been used in theories of elementary particles, condensed 
matter, accelerator and plasma physics, oceanographic 
and atmospheric sciences, and optics [16 . Central to our 
work is the recognition that any linear passive symplec- 
tic transformation S acting on the canonical Hermitian 
operators x^y^Px^ and Py (whose only nonvanishing com- 
mutators are [x,j9a;] = [y^Py] = ^A), is associated with a 
unitary operator U (S) generated by the group 
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Here, \ = 1/k and wq is the characteristic width of the 
spatial modes Um^n in which the analyzed wave field is 
to be decomposed. The set ([T]) satisfies the usual SU(2) 
algebra >C^] = iSahcJ^c h,c = x, y, z), with A/" being 
the only commuting generator in the group and jCz de- 
scribing real spatial rotations on the transverse xy plane 
(it is proportional to the orbital angular momentum op- 
erator along the wave propagation direction [14 )• We 
can therefore represent the most general passive unitary 
operator U (S) associated with by a single exponential 
of linear combinations of any of the above generators fiE\ 
U{S) = exp[— z(^+ A/' + ■ C)]^ with real scalar and 
vector parameters. We show below how U{S) can be 
implemented using simple optical elements. 

The method proposed here relies on exploiting A/", to- 
gether with specific combinations of the generators Cx 
and jCy, to construct two commuting unitary operators 
from which the associated symplectic matrices and their 
experimental implementation can be found relatively eas- 
ily. To this end, let |m,n) denote the (pure) mode 
states of order m + n > (in position-representation 
(x,?/|m,n) = Um,n)- We first impose |m,n) to be eigen- 
states of both JV and the operator Co^^p = Ur ■ where 
= (cos (f sin 0, sin cp sin ^, cos 0) may be conceived as 
a radially oriented unit vector in the orbital Poincare 
sphere [171 HE]- The eigenstates |m,n) depend on the 
choice of and (p and fulfill A/'lm, n) = [{m + n)/2] |m, n) 
and Ce^ip\m^n) = [(m — n)/2] |m, n). For instance, the 
eigenvectors of operators Cx and Cz are the HG and 
LG modes, respectively [14[|18]. Measurements of lip) = 
J2m n ^m,n\^^ ^) will iuvolvc the actiou of the unitaries 

Uj\f — e~'^^+-^ and Ue^ip = e"*'^-^^ '^, upon variation 
of parameters and (/)_, which can be externally con- 
trolled. Notice that Uj\f is connected with the Gouy 
phase [19], whereas Uq^^ describes ^_-angle rotations 
about Ur, thereby changing the mode superpositions. 

In order to extract the complete spectrum Vm,n of a 
coherent electromagnetic scalar field, an optical scheme 
is further developed. We remark that for quantum me- 
chanical matter waves (e.g. Bose-Einstein condensates), 
a conceptually similar approach to analyze their spa- 
tial structure should currently be feasible by exploiting 
atom optics: interferometers |20j, beam splitters in atom 
chips [21], focusing and storage in resonators [2^ , and 
conical lenses [23]. We use here a Mach-Zehnder inter- 
ferometer with built-in refractive components performing 
and Ue^ip' sets of spherical and cylindrical thin lenses. 
For ease of operation, it is desirable to vary (j)± in a way 
that minimizes the displacements of the optical elements. 
We explicitly show that all the required transformations 
in the interferometer can be achieved solely by rotation 
and variation of the focal lengths of the lenses. Trans- 
lations of the lenses are not necessary, although it may 
turn out to be more convenient in certain instances to 
carry finite displacements in some of them. In any case. 



the arms of the interferometer always remain fixed. 

Detection of the light intensity difference AI = Ib — Ia 
at the two output ports (A and B) of the interferometer 
provides the data for reconstructing the weights Vm,n- 
One has AI oc mU^^Ul^Uc + UlU^UeM). where 
Uc represents the unitary operator for the compensating 
lens system in the complementary arm of the interfer- 
ometer. For clarity, let us assume that Uc comprises 
a similar lens set as the one for Uj\fUe,if (we will show 
later that this assumption can be removed), with equiv- 
alent parameters 0^ and The mode spectrum and 
the output intensity difference are directly connected by 
a double Fourier-like transform 
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where the proportionality constant equals (IGtt^)"^ if 
m = n = 0, and (Stt^)"^ otherwise. Note that AI pos- 
sess the symmetry properties A/((/)+ ± 27r, (j)- ± 27r) = 
A/(^+ ± 27r, (j)- ^ 27r) = A/(^+, ^_). Our scheme non- 
trivially generalizes that of Ref. [11 , aimed to reveal the 
quasi-intrinsic nature of the orbital angular momentum 
degree of freedom for scalar waves. There, by means of 
the measurement U = e"*"^^^, implemented with Dove 
prisms, the azimuthal index spectrum of spiral harmonic 
modes could readily be accessed. Here, we explore the 
entire transverse mode space expanded by the indices m 
and n. Hence, the combined effect of the two nonequiv- 
alent measurements and Ue,(^ is crucial, and would 
allow us to measure, for instance, either the full spectrum 
of HG or LG modes. 

To show an explicit application of the above approach, 
we proceed with the characterization and design of an 
optical system that enables the reconstruction of the HG 
mode spectrum Vn^.uy of a light beam {ux^riy > 0). In 
this scenario, one needs to consider the unitary operators 
and t/6>=7r/2,(/)=o = e"*"^-^^. The path from these 
unitaries to the pursued setups resorts to the Stone-von 
Neumann theorem [I3l [HI [15]; it yields the symplectic 
matrices 5+ and S- which contain all systems informa- 
tion 
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where c± = cos(0±/2), s± = sin(^±/2), and zq = 
Wq/{2X) is the Rayleigh range. From matrices S± one 
can then obtain the corresponding integral transforms 
that govern the propagation (along the z-direction) of 
any input paraxial wave ip{x^y) traversing each system. 
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FIG. 1: (Color online) (a) Optical scheme for measuring the 
HG mode spectrum. Transformations 5*+ and S- consist of 
symmetric sets with three fixed tunable-focus spherical lenses 
and three pairs of rotary cylindrical lenses, respectively. The 
compensating system involves two settings with four and two 
spherical lenses, respectively, (b) Detail of one of the assem- 
bled pairs of cylindrical lenses required for S-. (c) and (d) 
Operation curves for 5*+ and respectively (see text). 



These transforms lead, in our case, to the integral kernels 
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The output wave profiles after each S± result from 
^±{x^y) = J dx'dy' K±{x^y; x' ^y')ip{x' ^y'). Kernels Q 
exhibit the structure of those for the fractional Fourier 
transform [24] . Integral transforms associated to unitary 
operators generated by Cy (rather than Cx) have been 
formulated [25]. They make possible the full meridian- 
rotation on the orbital Poincare sphere and, in particular, 
the reciprocal conversion between HG and LG modes. 

The two proposed optical systems for S± are symmet- 
ric and illustrated in Fig.jlja). For 6*+, three fixed spher- 
ical lenses with varying radii of curvature Rj (linked to 
the focal lengths fj by Rj = {fij — l)/j, nj being the 
refractive indexes), and placed at equal distances zq^ are 
sufficient. They fulfill Ri = R3 = zo/[l - cot((/)+/4)] 
and R2 = zo/(2 — 5+). Tunable-focus liquid crystal 
spherical lenses controlled by externally applied voltages 
have been demonstrated, displaying a wide range of focal 
lengths |26J . Figure [TJc) shows the operation curve de- 
scribing the dependence between Ri , R^ and R2 needed 
to cover the interval tt < 0+ < Stt. These values can be 
attained with the lenses of Ref. [2^ . For S- , three pairs 
of cylindrical lenses are required [see Fig. [ija)]. Each 
pair of assembled cylindrical lenses is rotated in a scissor 
fashion with a = — /3 [Fig. [TJb)]. The corresponding an- 
gles satisfy ai = as = (tt — l])/4 and a2 = (Stt — ^_)/4. 



The operation curve in Fig. [TJd) represents the variation 
of the rotation angles in accordance with the constraint 
imposed by cot(0_/4) = — 2sin(r^/2). This equation is 
satisfied when tt <</>_< Stt and — (tt/S) < Q < (tt/S), 
which lead to < 0^2 < (7r/2) and (tt/G) < ai < (tt/S). 
The radii of curvature of the cylindrical lenses are Ri = 
R2 = R5 = Re = zo/2 and R3 = R4 = ^o/4. The 
distance between consecutive pairs is 2^0/2. With this 
scheme, the covered values for G [7r,37r]. The fact 
that both and (/)_ are restricted to the interval [tt, Stt], 
rather than to [0,47r], does not constitute a fundamental 
limitation. It can be circumvented by employing the com- 
pensating system [see Fig. [ija)] and the symmetry prop- 
erties of A/. If two sequences of measurements are made, 
each having a different compensating system performing 
transformations with <p'_^ = (pL = (identity matrix) and 
0+ =0,0^ = 27T (minus identity matrix), respectively, 
then Eq. ([2| can be cast (m and n ^ Uy) as 

P„^,n, cx / / AJ(0+, - 27r) 

J TT J TT 

'{fix + ny)(l)^ + (n 
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The proportionality constant is (Stt^)"^ if Ux = riy = 0, 
and (47r^)~^ otherwise. The integration intervals in ^ 
now display the accessible ranges for (j)±. The first 
and second sequences of measurements can be carried 
out with a compensating system made of four and two 
(all identical) spherical lenses, respectively. By properly 
choosing their focal lengths, it is not necessary to displace 
the S± systems nor the interferometer arms. Moreover, 
measuring the LG spectra would only require embedding 
the system S- between two fixed, mutually orthogonal, 
cylindrical lenses, with /S+ remaining invariant. 

To demonstrate that our data analysis scheme is fea- 
sible and does not demand a large number of measure- 
ments for each (j)±^ we have numerically simulated the 
transformation and processing of several input coherent 
waves. Figure [2] depicts the light profiles entering into the 
interferometer: strongly astigmatic Gaussian [Fig. [2|a)], 
hexapole necklace [Fig. |2|d)] and anisotropic multiring 
[Fig.|2|g)] beams. Their exact HG distributions together 
with those retrieved from Eq. ^ are also displayed. To 
simulate the evolution of the various beams, the ker- 
nels Q were used to calculate AI in Eq. ([5| for 10 dif- 
ferent values (in equal increments) per each ^± G [tt, Stt]. 
The reconstructed and exact weights agree well (compare 
right and central columns in Fig.[2|. In analogy with the 
Whittaker- Shannon sampling theorem in Fourier analy- 
sis 0, input beams that are mode-band-limited can be 
exactly reconstructed via our scheme [compare Figs. [2|e) 
and|2|f)]. That is, if the sampling increments for are 
smaller than the inverse of the highest contributing mode 
numbers, reconstruction will be exact. Misalignment of 
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FIG. 2: Mode spectrum analysis. Exact and reconstructed 
[via Eq. ^] HG mode weights (central and right columns, 
respectively) from the input beam profiles (left column). 



the optical elements is expected to be the main source of 
errors. The effect on the weights Vn^.Uy-, due to lens dis- 
placements (tolerances) 5 with respect to the beam axis, 
introduces a correction term ~ ^ which is smaller 

than 1% for typical values 5 < 10/im and wq > 100/im. 

In conclusion, we have proposed a spatial spectrum an- 
alyzer for waves with arbitrary profiles based on a simple 
setup of lenses in a Mach-Zehnder interferometer. Mea- 
suring the mode spectrum of optical beams is fundamen- 
tal not only to achieve control of laser beam profiles but 
also to access information encoded in this degree of free- 
dom of the electromagnetic field, in the context of spatial 
multiplexing. The generality of the presented theoretical 
analysis encompasses different mode bases and suggests 
the measurement strategy also for matter waves. The 
specific set-up for HG spectra is presented with details 
in view of a more immediate experimental realization. 
Moreover, our scheme is compatible with mode analy- 
sis at the single-photon level, and could be of great use 
for quantum information applications. Retrieval of the 
complete spatial mode spectrum of any monochromatic 
optical field would become an attractive tool for signal 
analysis and processing when combined with novel holo- 
graphic recording materials designed for modal tailoring 
and multiplexors. In this respect, amorphous photopoly- 
merizable glasses, exhibiting high refractive index modu- 
lation in experimental demonstrations of the optical Pen- 
dellosung effect [27J, are well suited for this purpose. 
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